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ABSTRACT 

It is now well established that changes in the X-ray spectral state of black hole low- 
mass X-ray binaries are correlated with changes in the radio properties of those sys- 
tems. Assuming radio power is a proxy for jet power, we can say that the jet is 
continuously present in the hard state and undetectable (and therefore weaker) in the 
soft state. Since the different accretion states are also generally assumed to be associ- 
ated with different disc geometries - the hard state with a hot, thick flow, and the soft 
state with a cold, thin disc - we investigate the possibility that these two phenomena 
are linked; i.e., that the difference in disc geometry is the cause of the difference in 
observed jet power. We do this by comparing various measures of jet power in nu- 
merical simulations of accretion discs of differing temperatures and thicknesses. We 
perform these simulations using the general relativistic magnetohydrodynamic code 
COSMOS++ and a newly added cooling function, which allows us to regulate the disc 
scale height H/r at different radii. We find no apparent correlation between the disc 
scale height and jet power whenever we normalize the latter by the mass accretion 
history of each simulation. We attribute this result to the role that the "corona" plays 
in confining and accelerating the jet (our corona may also be considered a failed MHD 
"wind"). The properties of the corona do not vary significantly from one simulation 
to another, even though the scale heights of the discs vary by up to a factor of four. If 
this holds true in nature, then it suggests that the correlation between spectral state 
and jet power must be attributable to some other property, possibly the topology of 
the magnetic field. Alternatively, it could be that the corona disappears altogether in 
the soft state, which would be consistent with observations, but has so far not been 
seen in simulations. 
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1 INTRODUCTION 

Black hole low-mass X-ray binaries (LMXBs) exhibit ra- 
dio emission (associated with jets) whose properties change 
with the observed X-ray spectral and, to a less well de- 
termined extent, tempor al properties of the accretion disc 
l|Fender et al.| [2004. 2009). Multiple systems have now been 
observed to follow similar evolutionary tracks in X-ray 
hardness-intensity diagrams, with their radio properties 
correlating with the position of the so urce in this dia- 
;ram: e.g. GX 339- 4 (|Belloni et all 120051). XTE J1650-50 
Corbel et al.1 120041 1 . and Cygnus X-l (|Wilms et al.ll2006l ). 
The salient details are as follows: Whenever the system is 
observed to be "hard," there is a steady, compact, optically 
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thick radio jet. As the X-ray luminosity increases during an 
outburst, the radio luminosity likewise increases. At some 
point, the X-ray luminosity stops increasing and the spec- 
trum begins to soften. The spectrum continues to soften 
until the hard component effectively disappears and a pure 
"soft" state is reached (a weak power-law component may 
still be observed in the soft state). Radio emissions consis- 
tent with a jet have so far not been observed in this state, 
although the data cannot ru le out an intermittent, optically 
thin, jet ll Fender et al.|[2009h . In one of the most extreme ex- 
amples, deep observations of the persistent soft-state source, 
4U 1957+11, suggest the soft-state radio luminosity must be 
at least 300 times weaker t han the expected ha rd-state lumi- 
nosity for the same source (|Russell et al.|[201ll ) . After reach- 
ing the soft state, the system generally fades while remaining 
soft. As it fades, the disc component scales as I a T 4 , indi- 
cating a constant emitter size. Eventually the system returns 
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to the hard state at a significantly lower X-ray luminosity 
than the initial hard-to-soft transition. 

In this paper we focus on the jet dichotomy between the 
hard and soft states. In order to structure our argument, we 
take the view that almost any model of jets from black hole 
accretion discs can be boiled down to four main "ingredi- 
ents:" 1) an energy reservoir that can be tapped, likely the 
rotational energy of the black hole (although possibly the 
disc); 2) a magnetic field configuration capable of tapping 
that energy, which almost certainly means the black hole or 
inner disc must be threaded by a net poloidal magnetic flux; 
3) a confining medium to help accelerate the jet and convert 
the magnetic energy carried in the Poynting flux into kinetic 
energy of the plasma; and 4) some mechanism to introduce 
mass into the jet since a purely magnetic jet could not ra- 
diate. Note that the first two ingredients are necessary and 
sufficient for a jet (or at least an outflow) to exist, while 
the latter two are important for it to be observable as a jet. 
Having recapitulated jet models in this way, we ask which of 
the ingredients are likely to change dramatically when the 
disc switches from one spectral state to the other. 

Since the changes to jet power we are concerned 
with are found within individual sources, they can only 
be attributable to properties that can change on time- 
scales of a few days, since that is the observed time- 
scale of state transitions (see, for example, the review by 
iMcClintock fe Remillardl 2006). This eliminates black hole 
rotational energy (ingredient #1) from consideration, since 
the mass and angular momentum of the black hole are es- 
sentially constant over such time-scales. The next three in- 
gredients all remain plau sible, however. 

Nu merical work by iBeckwith et al.l (|2008l ) [and more 
recently iMcKinnev et al. ( l|2012D ] shows quite nicely that a 
change in magnetic field topology (ingredient #2) from a 
poloidal to a toroidal configuration leads to dramatically 
weaker or transient jets. This is because, with a toroidal 
field, the black hole is not fed a consistent poloidal flux, 
which it needs to build and maintain a strong jet. There- 
fore, a global field change could explain the observed jet 
dichotomy if the disc always maintains a net poloidal flux 
around the black hole in the hard state, but not in the 
soft state. An argum ent along these lines was made by 
llgumenshchevl (2009). However, this does not account for 
why the disc itself undergoes a significant change during 
the state transition, as in most respects in numerical sim- 
ulations, the properties of the disc are the same regardless 
of magnetic field topology (magnetically driven turbulence 
within discs is larg ely insensitive to the net magnetic flux, 
unlike jets) (|Beckwith et al.ll200Sf ). This implies that some- 
thing else besides the field topology must also change during 
the state transition in order to explain the different appear- 
ance of the disc. In our opinion, this need for two seemingly 
unrelated changes to happen nearly simultaneously every 
time a state transition occurs presents a significant hurdle 
for this picture. Therefore, for now, we also eliminate ingre- 
dient #2 from consideration. 

That leaves ingredients #3 and #4 as candidates. First 
we point out that these two ingredients could be linked: 
the same medium that provides confinement for the jet 
might also be responsible for mass loading it. This would 
be especially true if there were some geometrically thick, 
centrifugally-supported structure surrounding the jet, per- 



haps som ething like an advect ion-dominated accretion flow 
fADAF) dNaravan fc Yil 119941 ). The truncated disc model 
(|Esin et al.lll997l ) argues that just such a structure should 
be present in the hard state. Furthermore, according to 
this model, this structure is supposed to disappear in the 
soft state, to be replaced by a geometri cally thin Shakura- 
Sunyaev disc IjShakura fc Sunvaevll 19731 ). which extends all 
the way down to the innermost stable circular orbit (ISCO). 

In this paper, we explore whether this same geometri- 
cal argument can explain why jets are observed in the hard 
state but not the soft. Our basic premise is that the jet (in 
terms of ingredients #1 and #2) is actually present in both 
accretion states, but so much less luminous in the soft state 
(because ingredients #3 and #4 are missing) that it has yet 
to be observed. It may even be that the "jet" (fast outflow 
with a narrow opening angle) converts into a "wind" (slow 
outflow with a wide opening angle) when the system transi- 
tions from the hard to the soft state. Such a picture would 
be consistent with the recent discovery in LMXBs of wide- 
angle winds that are ob served in the soft state but not the 
hard (|Ponti et al.ll2012T ). 

The recent of incorporation of a new cooling function 
into the computational astrophysics code COSMOSH — h has 
presented us with an opportunity to test our idea numeri- 
cally. We perform a series of general relativistic magnetohy- 
drodynamic (GRMHD) simulations of accretion discs that 
differ only in their scale heights. We then measure the out- 
put jet power and look for a correlation between this and 
the disc thickness that would be analogous to the observed 
correlation in nature. Since we carefully track the mass ac- 
cretion history of each simulation and treat the backgrounds 
identically, we can be confident that any differences we ob- 
serve are due solely to the geometry of each disc (or equiv- 
alently the disc temperature, which controls the geometry). 
Numerical limitations in treating the mass loading of the 
jet prevent us from being completely quantitative with our 
results; nevertheless, we can robustly test our basic idea. 

Since this is the first GRMHD study of jets from rel- 
atively thin discs (H/r < 0.1), the results should also be 
of some general interest. Along with testing our premise 
about the jet dichotomy in LMXBs, this work also extends 
our understanding of jets from black hole accretion simula- 
tions. Previous studie s of jets from such simu l ations (e.g 
De Villiers et all 120051 : lHawlev fc Krolikl 120061 ; IMcKinnev! 
20061 ). did not include cooling functions. As a result, the 
scale height was rather large (H/r > 0.1) and could not be 
varied independently of other parameters. The only other 
study we know of that teste d jet properties specifica l ly as a 
function of disc thickness is iMoscibrodzka fc Progal (2009), 
although their discs were embedded inside an in-falling, 
quasi-spherical envelope and they controlled the disc thick- 
ness by adjusting the adiabatic index F of their gas. Our 
work is distinct from theirs in both respects. 

Our paper is organized as follows: In [J2]we describe the 
numerical scheme used in this work, including how the sim- 
ulations were initialized and the diagnostics used. In [J3]we 
present our main results, comparing jet power in simulations 
with different disc scale heights. We find that, contrary to 
our arguments, there is very little correlation between the 
output jet power and the disc scale height when one normal- 
izes for the mass accretion history. In [|4]we speculate why 
our numerical results did not support our hypothesis. Fi- 
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nally, in Sj5]we present some concluding thoughts. Through- 
out this paper we use units such that GM = c = 1; thus, 
units of length and time are tg = GM/c 2 and to = GM/c 3 , 
respectively. 



2 NUMERICAL METHOD 

We solve the following set of conservative, general relativistic 
MHD equations, including radiative cooling, 



dtD + d^DV 1 ) = , 

d t £ + di(-V=g n) 
d t s j + d i (^Ti) = 



-g n r t t + Aut 



-g Au-j 



d t B 3 +d,(B J V z -B Z V J ) = 



(1) 
(2) 
(3) 
(4) 



where D = W p is the generalized fluid density, W = ^J—gu 
is the generalized boost factor, V 1 — u l /u is the transport 
velocity, u M = g^ v u v is the fluid 4-velocity, g^ v is the 4- 
metric, g is the 4-metric determinant, 

£ = -V=gTt 

= -Wut{ P h + 2P B )-v^g-(P + P B ) + V^gb t b t (5) 

is the total energy density, h — 1 + e + P/ p is the specific 
enthalpy, e is the specific internal energy, P is the fluid pres- 
sure, Pb is the magnetic pressure, T£ is the stress-energy 
tensor, A is the cooling function (described in detail in i )2.3|) . 
and 



-gTJ = Wuj(ph + 2P B ) 



(6) 



is the covariant momentum density. With indices, V indi- 
cates the geometric connection coefficients of the metric; 
without indices, V is the adiabatic index. For this work we 
use an ideal gas equation of state (EOS), 



P=(r-l)pe 



(7) 



with T = 5/3. 

There are multiple representations of the magnetic field 
in our equations: 6 M is the magnetic field measured by an 
observer comoving with the fluid, which can be defined in 
terms of the dual of the Faraday tensor 6 M = u v * F M " , and 
= y/~gB :> is the boosted magnetic field 3-vector. The 
un-boosted magnetic field 3-vector B % = *_F Ml is related to 
the comoving field by 

B l = u l b l - u'b 1 . (8) 

The magnetic pressure is Pb = b 2 /2 — fo M 6 M /2. 

Equations JT]) - ((4} are solved using the new High 
Resolution Shock Capturing (HRSC) scheme recently 
added to the COSM OSH — h computational astrophysics code 
l|Fragile et alJl2012l ). 



(^centre = 45rG), and the power-law exponent (g = 1.54) 
used in defining the specific angular momentum distribution 
(and thus influencing the initial thickness), 



g±<i, + tgu 



9/2-1 



(9) 



t-9<H> + t- 2 9t4> 

We follow the procedure in IChakrabartH ljl985h to solve for 
the initial, hydrodynamic-equilibrium state of the torus. 

To this torus we add weak poloidal magnetic field loops 
in order to seed the ma gneto-rotational instability (MR1) 
( Balb us fc Hawlevl[l99lh . The non-zero spatial components 
are B r = —daA v and B® = d r A v , where 



C(p - p cut ) 2 sin [2N log {r/S)] for p > 
for p < p cut 



pcut 



(10) 



N = 1 is the number of field loop centres, and 5* = l.lrj n is 
a scaling parameter. The parameter p cut — 0.5p ma x,o, where 
Pmax,o is the initial density maximum within the torus, is 
used to keep the field a suitable distance inside the surface 
of the torus initially. Using the constant C, the field is nor- 
malized such that initially /3 ma g = P/Pb /3ma g ,o = 10 
throughout the torus. This is to ensure that the initial mag- 
netic field is weak. 

In the background region not specified by the torus so- 
lution, we set up a low density non-magnetic gas, with pro- 
files p = 10~ 4 p max ,or -2 ' 7 , e = pe = 10" 6 p max ,o?"" 2,7 , and 
1 — ( r G/r) 4 ]. Numerical floors of the form 
and e noor = 10 



V r 

Pflo. 



= 10 7 pmax,0r 2 ' 7 



9 —2.7 

Pmax.or are 



placed on the density and internal energy density. This den- 
sity floor is sufficiently low and falls off fast enough radially 
that it is very rarely applied in the simulations, even inside 
the relatively evacuated jets. Instead, the most commonly 
applied floor, and the source of matter inside the jets, is 
the floor applied on the ratio (p + pe)/Pe. Whenever this 
quantity drops below 0.01, both p and pe are rescaled by a 
factor appropriate to maintain this ratio. This procedure is 
necessary as the code can not handle more extreme magne- 
tization. 

Our simulations are run for 7.5 orbital periods, as cal- 
culated at the pressure maximum, r ccn tr C . This corresponds 
to ~ 14200£g, as the orbital period is ~ 1900£g- We con- 
sider three different target thicknesses (H*/r = 0.04, 0.08, 
and 0.16), plus two models where the disc has both a thick 
(H*/r w 0.16) and thin (H*/r ~ 0.04) component with a 
specified transition radius r tr . The names and parameters of 
each of the simulations are presented in Table [2~T1 Note that 
simulation hl6rAll has a target scale height {H,/r = 0.16) 
that is thicker than the disc can achieve over the course of 
the simulation (the starting thickness is H/r w 0.08). This 
was intentional, as we wanted to have one "no-cool" refer- 
ence simulation. It is important, though, that the non-disc 
material is treated consistently in all our simulations. 



2.1 Initialization 

All of our simulations start with identical initial condi- 
tions, consisting of a gas torus orbiting within a low den- 
sity background around a Kerr black hole. The free pa- 
rameters that describe the simulations are the black-hole 
spin (a* = a/M = cJ/GM 2 = 0.93), where J is the an- 
gular momentum of the black-hole, the inner radius of the 
torus (rj n = 15rG), the radius of the pressure maximum 



2.2 Grid 

Our simulations are carried out in the modified Kerr-Schild 
coordinates (t, x\, x%, <f>) in 2.5 spatial dimensions (all 
three spatial components of vector quantities are evolved, 
although symmetry is assumed in the azimuthal direction). 
The grid is uniformly spaced in coordinates xi and xi, with 
1024 and 512 zones, respectively. All curvature, both real 
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Table 1. Simulation parameters 



Simulation 


H t /r 


n r /rG 


H/r^_ 


H/r^ 


'£. 


h04rAll 


0.04 




0.09 ±0.02 


0.07 ±0.03 


0.39 ±0.11 


h04rl5 


0.16 


15 


0.10 ±0.02 


0.09 ±0.02 


0.47 ±0.08 




0.04 






0.04 ±0.01 




h04r30 


0.16 


30 


0.10 ±0.02 


0.11 ±0.03 


0.48 ± 0.08 




0.04 






0.04 ±0.01 




h08rAll 


0.08 




0.06 ±0.01 


0.06 ±0.01 


0.60 ±0.12 


hl6rAU 


0.16 




0.11 ±0.02 


0.12 ±0.02 


0.52 ± 0.14 



a Average of profiles in Fig. [T] from ngco = 2.10tq to IOtq, the 
approximate light-cylinder radius. Errors represent one standard deviation 
(Iff) about this mean. 

^ Average of profiles in Fig.[l]from ngco = 2.10rG to &0tq or nsco to rtr 
and rt r to 40tq. Errors represent lc about this mean. 

c Averaged over the final 2.5 orbits of the simulations, 9500 ^ t/tc ^ 14200. 
Errors represent lc about this mean. 



and coordinate, is handled via the metric terms in Cos- 
MOSH — h. The key then is the transformations between the 
grid coordinates xi and X2 and the true spatial coordinates, 
in this case the Kerr-Schild r and 6 coordinates. These trans- 
formations are modeled after iMcKinnevI l|2006l ). The radial 
coordinate is given by 



r = Ri + Roe x 



(11) 



where Ri controls the radial concentration of zones towards 
tbh and n r affects the radial distance covered by each zone. 
We choose Ri = —3, Ro = tbh, and n r — 10, allowing 
our grid to cover a large enough radial extent to study the 
disc at small scales and jet at large distances with a reason- 
able number of zones. Because we use horizon-penetrating 
Kerr-Schild coordinates, we are able to place the inner ra- 
dial boundary of the grid inside the black hole horizon, 
tbh = + yl — of) = 1.37rc, thus isolating it from the 
physical domain of interest. The radial range extends from 
0.9 tbh to 10 re- We use outflow boundary conditions at 

both the inner and outer radial boundaries. 

The polar angle transformation is given by ( Mc KinnevI 

l20oeh 



■■X2 + -[l-h{r)} sin(2a;2), 



(12) 



where h(r) alters the angular concentration of grid zones, 
with values between and 1 concentrating zones near the 
equatorial plane and values between 1 and 2 concentrating 
zones near the poles. We use the expression 



h(r) = 2 - (2 - i 



- 



(13) 



\ r 0j 

with Qj = 0.7 and roj = 3. This concentrates grid zones 
near the equator at small radii and near the poles at larger 
radii. The angular index, 



<-h 



1 1 

— I arctan 

2 7T 



(14) 



controls the transition between the inner and outer radial 
regions. We choose rij = 20, rij = 120, and rij = 0.3. The 
polar angle covers the full range ^ 9 ^ tv, and we use 
reflecting boundary conditions along the poles. 



The advantage of this grid construction is that it allows 
us to resolve both the disc and the jet. Even for our thinnest 
discs, our resolution is sufficient to ensure that the fastest 
growing poloidal field MRI modes are well resolved [Amri = 
8Az, where va is the Alfven speed]. We also 
confirm that q mag = —g^^b^jPB ~ 0.3 — 0.4, as expected 
(Hawl ev et alj|201ll ). At large radii, because the grid zoning 
becomes nearly cylindrical around the pole, we are able to 
maintain a nearly constant resolution across the jet core, 
which also becomes nearly cylindrical at large distances. We 
are also able to maintain a reasonable number of zones across 
the interface between the jet and the surrounding material, 
which can be the corona, failed MHD wind, or "background" 
gas. 



2.3 Cooling Function 

Since the goal of this work is to study how jet power depends 
on disc thickness in GRMHD simulations, we need a way to 
control this par ameter. We u se a m odified form of the cooling 
function from iNoble et all (120091 ). which adds a radiative 
loss term to equations ((2} and ©. The version we use has 
the form 



° 

A = fie I Y - 1 



Y < 1 
1 < Y < 2 

Y > 2 



(15) 



where Q(r) = M 1/2 /(r 3/2 +aM 1/2 ) is the relativistic orbital 
frequency and Y — (T — l)e/T* is the ratio of the actual to 
the target temperature. The target temperature is estimated 
from the target scale height as 



H* 



(16) 



The key difference b e tween our cooling function and the 
one from lNoble et all (|2009h is that, in our form, the cool- 
ing function has a maximum value of A ^ fie, whereas in 
INoble et all (|2009h . the cooling continues to be proportional 
to some power of Y, meaning the cooling rate can be quite 
high for gas that has a temperature far from the target value. 
This is an important distinction; in the case where A tx Y q , 
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the gas in the jets cools quite strongly, and that cooling 
would vary from one simulation to the next as we change the 
target scale height. Since our goal is to compare jet power, 
these artificial variations are unacceptable. In our version, 
the background gas and jets are largely unaffected by the 
cooling function, and only the disc is significantly altered. 

As one final modification to the cooling function, we 
allow the target scale height to vary with radius. In this 
work we consider either a single target scale height or two 
target scale heights, one that applies at small radii (and 
allows the disc to be thick) and one that applies at large radii 
(and forces the disc to be thin). The purpose of having two 
target scale heights is to mimic the truncated disc model, 
where a thin accretion disc at large scales is truncated at 
some radius, inside of which the flow is significantly hotter 
and thicker. Note that, in either case, the cooling of non-disc 
gas is the same, so our jets should not be affected directly 
by this difference. 

Since it is conceivable that the cooling time of the gas, 
tcooi = e/A, could sometimes be shorter than the MHD 
timestep, Atymu ~ Ax/V, where Ax and V are the char- 
acteristic zone length and velocity, respectively, we allow 
the cooling routine to operate on its own shorter timestep 
(subcycle), if necessary. To prevent runaway loops inside the 
cooling function, we limit the cooling subcycle to 5 itera- 
tions, with a maximum change to the specific internal energy 
in each subcycle of Ae = O.le. 

As a final note on the cooling function, we emphasize 
that it is not intended to capture any particular physical 
process. It simply acts as a heat sink mechanism, allowing 
the disc to maintain a specified thickness. Thus, it gives us 
very good control over the main variable we want to test. 
Another nice feature is that, because there is no physical 
scale associated with this cooling function, our results can 
apply to black holes of any mass. 



2.4 Diagnostics 

To measure the scale heights of our discs, we use the follow- 
ing density-weighted expression 



H 

r 



I- 



-gp 2 (8~n/2f d8dcj}\ 1/2 



I- 



(17) 



Using this diagnostic, we find that the discs in our simula- 
tions generally match their target scale heights quite well 
(as shown in Table |2~T1 and discussed in ij3.ip . 

To define our "jets" we used the criterion b 2 /p > 
1, emphas izing that these are magne ti cally-driven jets 
(see also iMcKinnev fe Gammiel |2004 |Pe Villiers et all 
l2005h. powered by the Bla ndford-Znajek (BZ) mechanism 
l|Blandford fc Znaieldll977r ). However, we wish to emphasize 
that, at this time, it is not clear exactly how our numerical 
jets might correspond to what obs ervers are able to mea- 
sure. iNaravan fc McClintockl (|201 lh argue that the ballistic 
jet observed in the steep power-law sta te is the one assoc i- 
ated wit h the BZ mechanism (but see iFender et al.ll2010l ). 
However, INaravan fc McClintockl (|201ll ) fail to explain why 
the BZ jet would only be seen in one accretion state. In 
our picture, the jet is expected to be present in all accretion 
states, since the fundamental ingredients (rotation plus open 
magnetic field lines) are argued to be present in all states. 



The question is then shifted to whether or not the BZ jet is 
observable in all states. 

We use a few different measures of jet power or efficiency 
in this work. The first i s the same measure of jet efficiency 
as iTchekhovskov et al.1 i|201ll ), namely rj = (M — E)/(M), 
where 



(18) 



M(t) = -j v-gpu 

is the mass accretion rate and 

E{t) = f V=gT[ 



(19) 



is the energy flux, both taken at the black hole event hori- 
zon tbh. Angle brackets indicate a time-averaged quantity. 
The negative sign in equation (|18[) indicates that the mass 
flux is positive when rest mass is transported into the black 
hole. Similarly, E is constructed such that a positive value 
indicates a net flux of energy into the black hole. 

A more global measure of jet power can be made by 
summing the energy carried in the electromagnetic and fluid 
components of the jet (b 2 / p > 1) 



E {t) = I y/=gTtdrd9d4> 



(20) 



where T t ' FL = phuu t + P and T t ' EM = 2P B uu t +P B - b f b t . 
We also track the energy flux E separately for the fluid and 
electromagnetic components along the length of the jet. This 
is useful for tracking the conversion of magnetic energy into 
kinetic. Finally we track the density-weighted Lorentz factor 



7 



\ s f^gp^f^gTtu 1 d6'd<j)' 
I ^gpd6'd(p' 



(21) 



where d0' and dcj>' are restricted to those parts of the radial 
shell that contain jet material (fc 2 // 9 > !)■ 

We also need to define a "corona" in our simulations. 
Without doing more detailed radiative modeling, it is diffi- 
cult to know for certain what component of our simulations 
might correspond to the corona observed in the hard-state of 
LMXBs. For the purpose of this paper, we use the working 
definition that the corona is comprised of all of the bound 
(— ut < 1) material surrounding the disc. Most of this ma- 
terial has /3 mag ^ 1. It was originally part of the torus, but 
gets lifted out of the disc by magneto-centrifugal forces. For 
these reasons it might be more appropriate to call our corona 
a failed MHD wind. 



3 RESULTS 
3.1 Scale Height 

Fig. Q] shows the scale height (H/r) as a function of ra- 
dius, averaged over the final 2.5 orbits of each simulation. 
Even with time averaging, there is some variability in the 
scale height; this mostly represents regions where the cooling 
function is unable to keep up with local heating (near the 
event horizon) or where there are gaps in the disc such that 
the scale height is effectively measuring the corona and not 
the disc (around 15rc in a few simulations). Particularly in 
simulation h04rAll, the disc breaks up into distinct clumps 
(rings) of matter, causing the scale height to appear large 
in the intervening gaps. Nevertheless, the geometry in each 
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Figure 1. Disc scale height H/r as a function of radius, averaged 
over the final 2.5 orbits, for each of our models. 



10" 



■fei io 



10 



10 



■tei 10 



10 



-i— r-i— i— i 

h04rAll 
h04rl5 
h04r30 



h08rAll 
hl6iA.ll 




3000 



6000 



9000 



12000 



t/t. 



Figure 3. Plot of the electromagnetic (-Eem) an d fluid (-Efl) 
energy fluxes, both in arbitrary units, at the black hole event 
horizon. 
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Figure 2. Plot of the magnetic pressure in the funnel (r = 6tq, 
~ 0), plus the total pressure in the corona (r = 6r<3, 9 = tt/3), 
and disk (r = 6rc 9 = 7r/2), all in arbitrary units. 
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Figure 4. Plot of the mass flux M and energy flux E, both 
in arbitrary units, at the black hole event horizon and the jet 
efficiency 77 = (M - E)/{M). 



simulation is largely as intended (see Table 12.11) , including 
the cases where the scale height changes from a thick to a 
thin component (simulations h04rl5 and h04r30). 



3.2 Magnetic Jet 

In th e Blandford-Znajek (BZ) model (jBlandford fc Znaiekl 
1977), the power of the jet is predicted to depend on the an- 
gular velocity of the footpoint of the magnetic field and the 
strength of the poloidal magnetic field component, Pbz oc 
SI 2 bp. The magnetic field is assumed to thread the event hori- 
zon, so the angular velocity is JIbh = a/2rBH- Since a (and 
therefore SIbh) is the same in all our simulations, we expect 
the Blandford-Znajek power of our jets should only vary if 
bp differs. Fig. [2] shows that, indeed, the magnetic pressure 
inside the funnel, which close to the pole is dominated by 
the bp component, does vary proportionally with H/r. This 
dependence is also reflected in the electromagnetic flux at 
the event horizon, Sem(^ = ^bh), as shown in Fig. [3] 



However, if we account for the fluid flux and normal- 
ize by M, i.e. we plot a jet efficiency, then the dependence 
on H/r disappears, certainly at late time, as shown in Fig. 
[4] (see also Table [2~Tj) ; all of our simulations settle to sta- 
tistically equivalent jet efficiencies of 77 w 0.5. The same 
independence is seen if we consider globally integrated mea- 
sures of jet energy normalized by the total accreted mass. 
For example, J?em/AM and J5fl/AM show no discernible 
trend with H/r in Fig. [5] 



3.3 Energy Conversion 

Having demonstrated that all of our jets contain the same 
power when normalized by the mass accretion history, we 
consider the question: How efficiently is this magnetic en- 
ergy flux converted into kinetic energy flux within the jet? 
Here we expect confinement to play a role, and it does, al- 
though not exactly in the way we anticipated. Fig. [6] shows 
that the jet is confined out to a radius r ~ 150rc, not by the 
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Figure 5. Plot of the total energy carried in the electromagnetic 
and fluid components of the jet, both in arbitrary units, normal- 
ized by the amount of mass accreted onto the black hole. 
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Figure 6. Contours of the magnetically-dominated jet boundary, 
b 2 1 p = 1 (solid lines), the corona, — ut < 1 (dashed lines), and the 
approximate disc surface, p = 0.01p max ,o (thin, dotted lines), for 
simulations h04rAll (black) and hl6rAll (purple). All quantities 
are time-averaged over the final 2 orbits. 



disc, but by what we label the corona. Fig. [6] also illustrates 
that the size, or extent, of the corona is not strongly depen- 
dent on the scale height of the disc, since it appears nearly 
identical for our two most extreme simulations (h04rAll and 
hl6rAll). This is consistent with our intention in setting up 
the cooling function ( H2.3P ; only the disc, not the corona, 
should be affected by changes to the target scale height. 
Fig. [7] then shows that the Lorentz factor indeed grows over 
the region of confinement, but the growth is modest, with 7 
only reaching values of 3-4 by r = 15f>G . 

The rest of the acceleration actually occurs after the 
confinement from the corona cuts off (r > 200rG in Fig. 0. 
The additional acceleration comes from the conversion of 
magnetic energy flux to kinetic, as can be seen in the bot- 
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Figure 7. Plot of the density-weighted shell average of the jet 
Lorentz factor 7 (top), the fraction of jet flux carried in the elec- 
tromagnetic component SemA^em + ^pl) (middle), and the 
fraction of jet flux carried in the fluid component -Efl/(Sem + 
Eel) (bottom). 



torn two panels of Fig. [7] This happens most efficiently once 
the matter passes t hrough the outer fast magnetosonic sur- 
face (lEichlerlll993l ; iBegelman fc Lilll994l ; lTomimatsulll994l ; 
iMcKi nncv 2006|). Nevertheless, the modest acceleration pro- 
vided by confinement is crucial, since it is only with this ac- 
celeration that the jet is able to achieve sufficient velocity to 
pass through the final critical surface and efficiently convert 
its energy to bulk kinetic motion. 

Ultimately, in the outermost parts of our simulation do- 
main (r > 5000rG), we end up with a near equipartition be- 
tween the magnetic and kinetic energy fluxes and a Lorentz 
factor of ~ 7, consistent with previous numer i cal studies 
llDe Villiers et all 120051 : lHawlev fc Krolikl 120061 : iMcKinnevI 
2006). Most importantly, we again see there is very little 
difference in the behavior of the jets; certainly the variation 
is much smaller than the factor of ~ 4 difference in scale 
height. 



4 COLD CORONA 

Our two main arguments for why thicker discs might lead to 
more powerful observed jets were that thicker discs would 
presumably lead to more mass loading and would also bet- 
ter confine their jets. However, the results in the previous 
section demonstrate two critical flaws in this argument: 1) 
confinement appears to be provided by the corona, not the 
disc; and 2) the bulk of the acceleration happens after the 
confinement ends. Nevertheless, this apparent contradiction 
of our original hypothesis may point us in a more promising 
direction: if it is the corona that provides the confinement, 
then perhaps the critical aspect of the hard-to-soft transition 
that could explain the fading of the jet is the disappearance 
of the corona. X-ray spectral observations already strongly 
suggest that while the hard state emission is dominated by 
a corona, the soft st ate spectrum shows little or no e vidence 
for this component ( Mc Clintock fc Rem illard 200^). If the 
corona indeed disappears in the soft state, then the only 
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question that remains is if this also explains the fading of 
the jet. 

Interestingly, the corona has never been observed to 
disappear in any GRMHD simulations that we are aware 
of, although there is some evidence for it being weaker in 
simulations that hav e no net poloidal flux (see Fig. 6 of 
iBeckwith et al.ll2008T l. In nature its disappearance may be 
related to some sort of thermal instability or change in ion- 
ization level, neither of which can be captured using the 
rudimentary treatments of gas thermodynamics currently 
implemented in global GRMHD disc simulations. Neverthe- 
less, we tried to create a smaller co rona by conducting a few 
simulations that used the original iNoble et all (|2009h cool- 
ing function, such that the non-disc gas cooled to a much 
greater extent than in the simulations presented here. Even 
so, we still found no significant differences between the coro- 
nae of these alternative simulations and the ones shown in 
Fig. [U] This, perhaps, should not be surprising since these 
coronae are generated by magneto-centrifugal, not thermal, 
forces. If the coronae in nature are similarly driven, then it 
may yet be that a change in magnetic field topology is re- 
quired to explain both the disappearance of the jet and the 



5 CONCLUSIONS 

The fundamental conclusion of this work is that there is 
no correlation between disc scale height and jet power in 
GRMHD simulations. Varying the disc scale height by al- 
most a factor of four (0.04 < H/r < 0.16) produces no 
meaningful differences in the measured jet efficiencies nor 
bulk Lorentz factors. 

Our main reason for suspecting that the disc scale 
height might affect jet power was that we presumed that 
the disc played an important role in providing confine- 
ment. Based on our simulations, it appears this is not the 
case; rather, it is the corona that provides the confine- 
ment. This find i ng co ntradicts the conjecture put forth in 
IBeckwith et alj (|2009T ) that it is the disk pressure that sets 
the level of (magnetic) pressure in the funnel. Interestingly, 
we also find that the corona is not sensitive to the level of 
cooling applied in our simulations, and so, its properties do 
not vary appreciably from one simulation to the next. 

In the Introduction, we suggested that jets ought to be 
present in all accretion states, since the necessary ingredients 
(rotation plus open magnetic field lines) seem likely to be 
present in all states. We would then attribute the differing 
jet observations in each state to the degree of confinement 
provided. In light of our results, we must modify this picture 
to either suggest that the different jet behavior between the 
hard and soft states must be due to a collapse of the corona 
or a change of field topology (scenarios that could not be 
explored using the methods of this paper). 

As always, there is concern as to whether or not our 
simulations have run long enough to reach meaningful equi- 
librium states. This is especially pertinent for our 2D sim- 
ulations, as these can never truly reach a steady state. The 
reason is that, after a period of initial vigorous growth, 
the MRI in 2D simulations steadily decays because the dy- 
namo action that normally sustains it requires access to non- 
axisymmetric modes that are obviously inaccessible. There 



are other pathologies related to the forced axisymmetry of 
2D simulations, which may explain, for instance, the gaps 
that are observed in some of our disks in Figure [T] 

Another issu e is that our simulatio n s, like many pre- 
vious ones (e.g. iDe Villiers et all 120051 ; lHawlev fc Krolikl 
120061 ; iMcKinnev fc Blandfordll2009l ). only supply a limited 
amount of magnetic flux to the black hole. This restricts 
the range of jet efficiencies that we can hope to observe. 
If, on the other hand, a very large amount of magnetic 
flux is available to the black hole, the black hole flux can 
saturate, with magnetic pressure reg ulating the accretion 
proces s, as demonstra t ed re ce ntly by Igumcns hchev et al.l 
(I2003T): llgumenshchev) (|200ST) : iTchekhovskov et all l|201lD ; 
IMcKinnev et al. (2012). In this "magnetically choked" state, 
the jet power does depend on disc scale height, as the thick- 
ness of the disc plays a role in regulating how much flu x 
reaches the black hole (|Tchekhovskov fc McKinnevll2012T ). 
Since our simulations only supply a limited amount of flux, 
we are not able to observe this effect. 

Disc scale height may also play a role in regulating 
jet power by setting how much of the black hole horizon 
participates in the BZ mechanism; only that portion where 
b 2 1 p > 1 can power the jet. In this case, it would only be 
the scale height very close to the event horizon that would 
matter. Figure [S] indicates that this event horizon "covering 
factor" is nearly identical for our two most extreme simula- 
tions. Thus, we do not expect to observe this effect, either. 

Clearly these concerns leave room for further explo- 
ration of the topic of how disc scale height might effect jet 
power. Our results also suggest that further study of the 
coronae in GRMHD simulations is warranted. 
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